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ABSTRACT 

Barred galaxies are known to possess magnetic fields that may affect the properties of bar substruc- 
tures such as dust lanes and nuclear rings. We use two-dimensional high-resolution magnetohydrody- 
namic (MHD) simulations to investigate the effects of magnetic fields on the formation and evolution 
of such substructures as well as on the mass inflow rates to the galaxy center. The gaseous medium is 
assumed to be infinitesimally-thin, isothermal, non-self-gravitating, and threaded by initially uniform, 
azimuthal magnetic fields. We find that there exists an outermost a; i -orbit relative to which gaseous 
responses to an imposed stellar bar potential are completely different between inside and outside. 
Inside this orbit, gas is shocked into dust lanes and infalls to form a nuclear ring. Magnetic fields 
are compressed in dust lanes, reducing their peak density. Magnetic stress removes further angular 
momentum of the gas at the shocks, temporarily causing the dust lanes to bend into an "L" shape 
and eventually leading to a smaller and more centrally distributed ring than in unmagnetized models. 
The mass infiow rates in magnetized models correspondingly become larger, by more than two orders 
of magnitude when the initial fields have an equipartition value with thermal energy, than in the 
unmagnetized counterparts. Outside the outermost xi-orbit, on the other hand, an MHD dynamo 
due to the combined action of the bar potential and background shear operates near the corotation 
and bar-end regions, efficiently amplifying magnetic fields. The amplified fields shape into trailing 
magnetic arms with strong fields and low density. The base of the magnetic arms has a thin layer in 
which magnetic fields with opposite polarity reconnect via a tearing-mode instability. This produces 
numerous magnetic islands with large density which propagate along the arms to turn the outer disk 
into a highly chaotic state. 

Subject headings: magnetohydrodynamics — galaxies: ISM — galaxies: kinematics and dynamics — 
galaxies: nuclei — galaxies: spiral — ISM: general — shock waves 



1. INTRODUCTION 

One of the characteristic features of barred galax- 
ies is the existence of gaseous substructures such 
as a pair of dust lanes, a nuclear ring, and nuclear 
spirals, visible in optical images and radio maps 
(e.g.. Pease 1917; Sandagc 1961; Sanders & Huntlc^ 
19761: iButa fc Combes! Il99a IMartini et all I2003ail: 



Prieto et all l2005t [Martinez- Valpuesta et all 12006 
Comeron et al.l l2009t iHsieh et al.l 120111) . Dust lanes 



located at the leading side of the bar are interpreted as 
shocks in the gas fiows, as evidenced by the enhanced 
radio emissio ns and sh arp velocity jumps across them 
(see, e.g.. iButa fc Combe s 1996. and references therein). 
The curvature of dust lanes depends on the bar strength 
and aspect ratio, while largely independent of its 
pattern speed and galaxy mass (jKnapen et al.l 120021 : 
iComeron et~alll2009f ). Nuclear rings are the regions of 
high gas density and usually populated with numerous 
H II regions, indicative of recent starburst activities 
(e.g. iButil 119861: IHeller fc Shlosmanl |T99l IMaoz eT^ 
[20qitlMazzuca et al.ll2008l: IHsieh et al.ll2011[ ). They have 
a size typically of 1 kpc, but appear to be larger and 
wider in galaxies wi th less compact mass distribution 
(jMazzuca et al.|[201lD . Nuclear spirals located inside the 
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rings are believed to be a channel for gas infiows to the 
galaxy center, potentially powering active galactic nuclei 
(AGN) (e.g., [Sh losman et al. 1990; Regan fc MulchaevI 
fl99llKnapen et al.ll2000l: Ivan de Ven fc Fathil[20M 

Since the gas flows shaping bar substructures are 
intrinsically nonlinear, numerical simulations have been 
a powerful tool to study their fo rmation, evolu tion, and 
associat ed dynamics (e.g., ISanders fc Pre ndergas^ 
JAthanassoulal Il992bl: [Piner et al.l 1199^ 



1971 



Englmaier fc Gerhardl Jl997i iPa tsis fc Athanassoulal 
20001: iRegan fc TeubenI 2003. 2004: Ma cieiewskil 12004 



Ann fc Tiiakui^'2005|: iThakur et al...2009iy . In particular, 



Athana ssoula (1992^) showed that a non-axisymmetric 



bar potential induces dust-lane shocks at downstream 
from the bar major axis and that the shocks be- 
come more straight under a stronger bar. Dust lanes 
tend to become shorter and located closer to the 
bar major axis when t he ga s sound speed is larger 
(lEnglmaier fc Gerhardl Il997t IPatsis fc Athanassoulal 
120001 ). The shocked gas at dust lanes loses angular 
momentum and flows inward to form a nuclear ring at 
the position where th e external gravity is balanced by 
the c e ntrifugal forces (iShlosnian et al.lll990l : iPiner et all 
119951: IRegan fc Teubenll2003D . Using a smooth e d par - 
ticle hydrodynam i cs m ethod, lAnn fc Thakuri (|2005D 
and iThakur et al.l (|2009( ) showed that the formation of 
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strong nuclear spirals supported by shocks requires a 
large sound speed and the presence of a central black 
hole (BH), although their models were unable to resolve 
weak spirals in the nuclea r regions due to a small 
number of particl es (see also iMacieiew ski' '2004*) . 

Very recently, iKim et all ( 20121 hereafter Paper I) 
ran high-resolution simulations using the grid-based 
CMHOG code on a cylindrical geometry. P aper I found 
that t he original CMHOG code used by iPiner et aLl 
(|1995l ) contained errors in the evaluation of the bar 
forces, comp romising their numerica l results as well as 
others (e.g., IRegan fc TeubenI 120031: IMacieiewskil I2004D 
who unknowingly adopted the original code with the er- 
rors. Using a corrected version of the code. Paper I 
found that when the gas sound speed Cg is small, nu- 
clear rings are narrow and decoupled from nuclear spi- 
rals, making the mass inflow rates M toward the galaxy 
center extremely small (~ 10""* M© yr~^). On the other 
hand, strong thermal perturbations in models with high 
Cs make the rings broad, leading to M ^ 10~^ M© yr~^. 
Paper I also found that the shape and strength of nuclear 
spirals depend sensitively on the gas sound speed and the 
BH mass Mbh such that they are leading if both Cg and 
Mbh are small, weak trailing if Cg is small and Mbh is 
large, and strong trailing if both Cg and Mbh are large: 
nuclear spirals are readily destroyed by the ring material 
on eccentric orbits when Cs is large and Mbh is small. 

While aforementioned work improved our understand- 
ing on gas dynamics associated with bar substructures, 
they are without one of the most important ingredients 
of the interstellar medium (ISM), name ly magnetic fields 
that are pervasive in disk galaxies (e.g.. lBeck et al.|[l996l : 
[Beck 2009). By observing barred galaxies in radio polar- 
ized emissions, iBeck et all ()2005D found that magnetic 
fields are strong in dust lanes and nuclear rings, sug- 
gesting that they are dyn amically important in shaping 
bar substructures (see also lBeck et al.lll999l ,l200?). Some 
barred galaxies such as NGC 1365, NGC 1097, NGC 1365 
are observed to possess "magnetic arms" characterized 
by stro nger fields and lower d ensity than surrounding 
regions (jBecket all [2002ll2005l ) . In addition, magnetic 
fields exert torque that is much larger than gravitational 
torque, providing an efficient means to tr ansport gas 
from a nuclear ring inward to feed an AGN (^Beck et al.l 
120051) . Therefore, the effect of magnetic fields on struc- 
ture and evolution of barred galaxies should not be ig- 
nored. 

There have been a number of numerical studies on 
magnetic field distributions in barred galaxies. Depend- 
ing on how magnetic fields are treated, they can be cat- 
egorized into two grou ps: (1) those based on meaii-field 
dyna i no theories (e.g., lOtmianowska-Mazur et al.|[T997l 
I200I IMoss et al.lll998l 119991 12001L I2007D and (21 those 
using full magnetohydrodynamic (MHD) simulations 
(iKulesza-Zvdzik et al.l [200I \20Wt iKulpa-Dvbet et al.l 
|2011|) . Although numerical models from the first group 
are successful in obtaining synthetic polarization maps 
and overall field morphologies comparable to observa- 
tions, they rely on parameterized turbulent terms in 
the induction equation that are uncertain. More impor- 
tantly, the mean-field dynamo models take velocity fields 
from hydrodynamic or iV-body sticky-particle simula- 
tions and evolve magnetic fields passively, without con- 



sidering the back reaction of magnetic fields to the gas. 
On the other hand, MHD models in the second group 
naturally handle gas responses to the embedded mag- 
netic fields. For instance, iKulesza-Zvdzik et al.l (|20100 
found that dynamical effects of magnetic fields lead to 
the formation of low-density magnetic arms in the regions 
outside the bar, which is clearly unseen in the mean- 
field dynamo models (see also lKulesza-Zvdzik et al.ll2009l : 
|Kulpa-Dybcl et al. 201l|). Even without the prescrip- 
tion for the mean-field dynamo, magnetic energy in these 
MHD models grows due to shear and compression by a 
factor of ab out 15 relative to the initia l value over a 1 
Gyr period (|Kulesza-Zvdzik et al.l[2009l ). 

Because the MHD models mentioned above focused 
on magnetic structures in the outer regions, however, 
the effects of magnetic fields on bar substructures and 
mass infiow rates have yet to be explored. Can mag- 
netic fields make dust lanes stronger? Do they make 
nuclear rings smaller or larger? Do nuclear spirals sur- 
vive in the presence of magnetic fields? How do the 
mass infiow rates change as the field strength varies? 
In order to address these questions, we in this pa- 
per run hig h-resolution MHD simulations using the 



per run nig n-resoiution MilJJ simulations using tne 
Athena code (iGardiner fc Stonel[20()5l: IStone et aLlMisI : 



IStone fc Gardineri I2009n ~ This work is a straightfor- 
ward extension of Paper I by including magnetic fields. 
Our models are two dimensional, assuming a razor-thin 
disk. This allows to better resolve the in-plane direc- 
tion than in the three-dimensional models considered by 
IKulesza-Zvdzik et al.l ()2010[ ). On the other hand, our 
two-dimensional models are unable to capture the dy- 
namical consequences of Parker instability and other pro- 
cesses that involve the vertical direction. In addition to 
studying the magnetic effects on substructures in the bar 
regions, we identify an MHD dynamo mechanism that 
amplifies magnetic fields, eventually leading to magnetic 
arms in the outer regions. 

This paper is organized as follows. In Section 2, we 
describe our numerical methods and model parameters, 
and present numerical resistivity of the Athena code. In 
Section 3, we make a comparison between the results of 
the Athena and CMHOG codes on hydrodynamic models 
with identical parameters. In Section 4, we present the 
results of MHD simulations on the bar regions, focusing 
on the effects of magnetic fields on bar substructures and 
mass inflow rates. In Section 5, we study evolution of gas 
and magnetic fields in the outer regions. In Section 6, we 
conclude with a summary and discussion of our numerical 
results. 



2. MODELS AND METHODS 

We consider a magnetized gaseous disk and study its 
responses to an imposed non-axisymmetric bar poten- 
tial. Gas dynamics without the effect of magnetic fields 
was presented in Paper I. As in Paper I, the disk is ini- 
tially uniform, isothermal, vertically-thin, and non-sclf- 
gravitating. The bar is assumed to rotate rigidly about 
the galaxy center with a fixed pattern speed Jib = i^bZ- 
We solve the dynamical equations in the frame corotat- 
ing with the bar in the z = plane. The equations of 
magnetohydrodynamics in the rotating frame integrated 
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u - V S = -SV • u, 



along the vertical direction are 

d_ 

dt 



(1) 



47rl] 



(V X B) X B 



f7gR-2J7b X u, 



dB 

'dt 



V X (u X B). 



(2) 



(3) 



Here, S, u, and Cg denote the surface density, veloc- 
ity in the rotating frame, and the sound speed in the 
gas, while B represents the midplane value of the three- 
dimensional magnetic field times the square root of the 
ratio of surface densit y to midplane volume density (e.g., 
IKim fc Ostrikeril2001[ ) . Note that the induction equation 
^ in the rotating frame has the same form as in the in- 
ertial frame (e.g., IKoldoba et al.l [20021 ). The velocity v 
in the inertial frame is obtained from v = u -I- RQ,b4'- We 
do not consider star formation and the associated gas 
recycling in the present work. 

The external gravitational potential, represented by 
$ext in equation and the associated rotation curve 
are taken identical to those in Paper I. Here we describe 
these briefly. The external potential consists of four com- 
ponents: a stellar disk modeled by a Kuzmin-Toomre 
disk, a central bulge modeled by a modified Hubble pro- 
file with cent ral dens i ty Pb ui, a non-axisymmetric bar 
modeled by a iFerrersI (|1887D spheroid with central den- 
sity pbar, and a central supermassive BH with mass Mbh- 
Without a BH, the rotation velocity Vc rises steeply with 
the galactocentric radius R ^1 kpc and attains a more- 
or-less constant value of ~ 200 km s~^ at i? ^ 6 kpc. 
The presence of a BH affects Vc near the central parts, 
resulting in Vc oc {Mbh/RY^^ at R<0.1 kpc. The bar 
has a patten speed of ilf, — 33 km s^^ kpc^^, which 
places the corotation resonance (CR) at Rqr — 6 kpc. 
In order to avoid abrupt responses of the gas due to a 
sudden introduction of the bar, we increase pbar linearly 
with time over one bar revolution time of 27r/rib = 186 
Myr, while keeping the net density Pbar + Pbui fixed. This 
ensures that shape of $oxt averaged along the azimuthal 
direction is not changed much with time. 

Our model disks initially have uniform surface den- 
sity So = 10 Mq pc~^ and isothermal sound speed of 
Cg = 5 km s~^. Magnetic fields are initially in plane, 
uniform, and purely azimuthal. The field strength is pa- 
rameterized by the dimensionless plasma parameter 
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where p and i?3D are the midplane values of volume den- 
sity and three-dimensional magnetic field strength. In 
terms of dimensional units 



B3B — 



4.6^G 

Vskms" 



10 Mq pc- 



1/2 



H 



100 pc 

(5) 



Table 1 

Model Parameters 



Model 


Mbh(Mq) 


/3o 


bh7MHD01 


4 X 10'' 


1 


bh7MHD03 


4 X 10'' 


3 


bhTMHDlO 


4 X 10'' 


10 


bh7HD 


4 X 10'' 


oo 


bhOMHDOl 





1 


bhOMHD03 





3 


bhOMHDlO 





10 


bhOHD 





oo 



-1/2 



where H = E/(2p) is the scale height of the gaseous disk. 
Note that the initial disk with constant /3o has a mass-to- 
flux ratio M/<^B = (7r/3oSo/8Cg)^/^i? that is increasing 
linearly with R. 

In order to focus on the effects of magnetic fleld 
strength and rotation curve at the central parts, we con- 
sider eight models with differing Mbh and /3o, as listed 
in Table [TJ The strength of magnetic fields is chosen to 
vary from sub-equipartition (/3o = 10) to equipartition 
(/3o — 1) values with thermal energy. We also consider 
hydrodynamic models with /3o = c» to make compar- 
isons with the results of Paper I. Models with a prefix 
bh7 possess a central BH with A/bh = 4 x 10'' Mq, 
with the f2 — k/2 curve attaining the local maximum 
and minimum at iJmax = 0.53 kpc and Rmin = 0.19 kpc, 
respectively (see Paper I). The bh7 models have only a 
single inner Lindblad resonance (ILR) at i?iLR ~ 2 kpc 
where Q. — k/2 — fi^. On the other hand, bhO mod- 
els with no BH have two ILRs at i?iiLR = 0.19 kpc 
and i?oiLR ~ 2 kpc, and the 17 — k/2 curve peaks at 
R^^^ = 0.53 kpc. We take Model bhTMHDOl with 
Mbh = 4 X 10"^ Mq and /3o = 1 as our standard model. 

We solve the time-dependent ideal MHD equations ((TJ- 
([3]) using a modified version of the Athena c o de in Carte- 
sian coordinates ([G ardiner fc Stond 120051 : iStone et"al] 
120081 : IStone fc Gardiner,, 2009^ . Athena utilizes a higher 
order Godunov scheme which conserves mass, momen- 
tum, and magnetic flux within machine precisions. It 
also provides several different schemes for integration in 
time, spatial reconstruction, and solution of the Riemann 
problem. We take the van Leer algorithm with piecewise 
linear reconstruction and first order flux correction. For 
a Riemann solver, we use the exact nonlinear solver for 
hydrodynamic models and an approximate HLLD-type 
nonlinear solver for MHD runs. 

Our simulation domain is a square box with side 
i = 30 kpc in each direction. We set up a uniform 
Cartesian grid with 4096 x 4096 zones over |?/| < 
L/2. The corresponding grid spacing is Ax — Ay — 
7.3 pc, about an order magnitude smaller than models in 
iKulesza-Zvdzik et al.l (|2010) . These high- resolution runs 
are necessary to explore properties of bar substructures 
in detail. We apply the outflow boundary conditions at 
the domain boundaries (i.e., at |a;| = L/2 or \y\ — L/2). 
To measure the mass inflow rates toward the galaxy cen- 
ter, we place a central hole with radius of 20 pc and deflne 
a density floor Sfioor inside the hole. At each time, we cal- 
culate the gas mass Min flown into the hole by integrating 
the excess interior density as Mjn = J^^^^^{Y,—J^iiooY)dxdy, 
'and then reset the interior density to Sfloor- The mass 
inflow rate is then obtained from M = Mjn/ At, where At 






X (kpc) 

Figure 1. Comparative snapshots of surface density in logarithmic scale from hydrodynamic runs with no BH using CMHOG (left) and 
Athena (right) at t = 300 Myr. The upper panels cover the corotation resonance, while the lower panels zoom in the central 1 kpc regions. 
The dotted ovals in the upper panels draw the outermost xi-orbit crossing the x- and y-axes at Xa = 3.6 kpc and yc = 4.7 kpc, respectively, 
while the short line segments indicate a slit perpendicular to the dust lane, along which density and velocity are measured in Figure [2k. b. 
The solid curves in the lower panels draw the X2-orbit with Xc = 1.2 kpc. 



is the computational timestep. In this paper, we present 
the results with Efioor = 10^^ M0 pc^^, but we checked 
that M is insensitive to the choice of Sfloor as long as it 
is taken sufhciently small. 

The inflowing gas may carry magnetic fields to the 
galaxy center, as well. In order to allow for B-field accre- 
tion, we let the magnetic fields vanish inside the central 
hole while preserving the divergence-free condition. This 
is achieved by calculating the vector potential at each 
time, cutting it off to a constant value inside the hole, and 
then recalculating the B vector from the modified vector 
potential. We note that taking B = Q arbitrarily inside 
the hole sometimes makes i3-fields bend abruptly near 
the hole boundary. In addition, density in the vicinity 
of the hole frequently becomes very small, greatly limit- 
ing the computational timestep by increasing the Alfven 
speed. In order to minimize these spurious effects, we 



utilize the non-ideal Ohmic dissipation term included in 
Athena by taking magnetic resistivity rj = 1 kpc^ Myr~^ 
in the circular region at i? < 40 pc, while the rest of the 
simulation domain is kept non-resistive. This results in a 
gradual dissipation of i?-fields in the region surrounding 
the hole, allowing us to run simulations for a long period 
of time. 

Finally, we remark on the numerical resistivity of the 
Athena code that we use. It is difficult to estimate the 
magnetic diffusivity of the real ISM accurately since it is 
very clumpy and diffusivity changes a lot between the dif- 
ferent phases with different temperature and different de- 
gree of ionization. If turbulence dominates the diffusion 
and dissipation of the magnetic fields, it is perhaps of the 
order of 7?T ~ IQ-^-lQ-^ kpc^ Myr"^ (e.g.,|Parker'1971|: 
iCamenzind fc LeschI 119941: iKulesza-Zvdzik et al.. ,200911 . 
The associated magnetic Reynolds number is Rm = 
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VcL/r]T ^ 10'' — 10^, much larger than unity, so that 
taking the Rm — >■ oo limit, as we do in the present work 
except for the small central regions, is a useful starting 
point. Nevertheless, we note that there is non-zero (but 
small) numerical diffusivity over the entire simulation do- 
main introduced by the numerical scheme. In Appendix, 
we use linear-amplitude magnetosonic waves to evaluate 
the coefficients of numerical viscosity as well as magnetic 
diffusivity of the Athena code. For Ax ~ 7.3 pc and the 
wavelength A ~ 0.5 kpc typical for the tearing-modes 
of magnetic reconnection occurring in the outer regions 
(see Section [5]), the numerical resistivity is estimated to 
be ??„ ~ 6 X 10""^ kpc^ Myr""', about three orders of 
magnitude smaller than ry-r- The associated magnetic 
diffusion time is Tmag = Ax^ jrin ~ 100 Myr, indicating 
that the non-ideal effect of numerical resistivity on the 
i?-field evolution can be important if the field strength 
changes considerably over the grid scale. 

3. COMPARISON BETWEEN CMHOG AND ATHENA 
RESULTS 

As mentioned above, Paper I explored hydrodynamical 
models of barred galaxies using the corrected version of 
the cylindrical CMHOG code tha t originally contai ned a 
bug in calculating the bar forces (jPiner et al.lll995D . On 
the other hand, the present calculations make use of the 
Athena code on a uniform Cartesian grid. Since the cur- 
rent models include two hydrodynamic models with the 
same parameters as in Paper I, it is interesting to com- 
pare the results from two different codes. This will check 
not only the reliability of the Athena runs but also the 
validity of the corrected bar forces in Paper I. Note that 
by employing a non-uniform grid, the CMHOG mod- 
els in Paper I has a higher (lower) spatial resolution at 
i? < 1.1 kpc (i? > 1.1 kpc) than the current Athena runs 
with a uniform grid spacing. 

Overall evolution of hydrodynamic models is described 
in Paper I: here we summarize the main features. As the 
bar potential is slowly introduced, initially-circular gas 
orbits are perturbed, forming a pair of overdense ridges 
at the downstream side of the bar major axis. As the 
bar potential grows further, the overdense ridges eventu- 
ally develop into shocks, appearing as dark dust lanes in 
optical images of barred galaxies. Gas loses a significant 
amount of angular momentum at the shocks and thus 
moves radially in, forming a nuclear ring at the position 
where the centrifugal force balances the external gravity. 
In addition, these shocks are curved and thus able to gen- 
erate vorticity. Vorticity grows secularly with time by the 
successive passages across shocks. When it achieves sub- 
stantial amplitudes, it causes the shock fronts to wiggle, 
producing small clumps with high vorticity. The whole 
process of clump formation is analogou s to the wiggle 
instab ility of spi ral shocks identified b y iWada h Kodal 
(|2004D (see also iKim fc Ostrikeii [2001) . These clumps 
are added to the nuclear ring. As clumps collide with 
each other, the ring becomes homogeneous gradually. 

Figure [T] compares the snapshots of Model bhOHD 
in the present work with those of Model csOSbhO from 
Paper I at i = 300 Myr when the bar substructures 
reach a quasi-steady state. The bar is oriented vertically 
along the j/-axis, and the gas inside (outside) the CR 
at i? = 6 kpc is rotating in the counterclockwise (clock- 
wise) direction. The images in the upper panels cover the 
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Figure 2. Comparisons between tlie Atliena and CMHOG results 
for tile distributions of (a) surface density and (b) velocity lij^ 
perpendicular and U|| parallel to the dust-lane shocks along the 
slit shown in Figure [T] (c) for the radial distribution of gas density 
averaged both azimuthally and temporally over t = 300 — 500 Myr, 
and (d) for the temporal variations of the mass inflow rates. In (c), 
Rxl = 4.7 kpc with an arrow marks the maximum radial extent of 
the outermost xi-orbit. 



whole region inside the corotation, while the lower panels 
enlarge the central 1 kpc regions. In the upper panels, 
dotted curves draw the outermost xi-orbit that cuts the 
X- and y-axes at Xc = 3.6 and yc — 4.7 kpc, respectively, 
and has a Jacobi energy of Ej = —1.24 x 10^ (km s~^)^ 
under our adopted external potential (e.g.. lAthanassoulal 
Il992al) . Inside this orbit, there are families of closed xi- 
and a;2-orbits that cross each other, resulting in shocks 
at the leading side of the bar major axis. As gas loses 
angular momentum at the shocks and moves toward the 
central parts, the bar regions becomes gradually evacu- 
ated. Outside the outermost xi-orbit, no closed orbit ex- 
ists: gas there simply moves on slightly perturbed, near- 
circular orbits, without inducing shocks. The gas moving 
in slowly due to the bar torque from outside is subse- 
quently trapped in the outermost -orbit, and stays in 
it. Therefore, the outermost cci-orbit acts as a barrier 
that inflowing gas cannot penetrate. Other than weak 
trailing spirals produced by the bar toque that emerge 
from the downstream side of the bar ends, the outer re- 
gions are almost featureless in hydrodynamic models. 

In the lower panels, solid curves draw an a;2-orbit with 
Xc = 1.2 kpc that approximately follows the nuclear ring. 
The morphological agreement between the CMHOG and 
Athena results is fairly good, although the ring in the 
CMHOG run is more clumpy and slightly larger. This 
is because the Athena run has a higher resolution at 
R > 1.1 kpc where dust lanes form. This makes the 
dust-lane shocks in the Athena runs better resolved, lead- 
ing to higher angular momentum loss and a smaller ring 
than in the CMHOG counterpart. In addition, the wig- 
gle instability of the dust lanes is more vigorous in the 
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Figure 3. Logarithm of the density distribution (left) and mai 
4 X lO'' M0 and /3o = 1 at t = 150 Myr. The shot line segment at 
dust lane, along which fluid variables are measured in Figure [4] 

Athena run since it grows faster at smaller scales. Hence, 
the Athena run produces more clumps that collide fre- 
quently, resulting in a smoother ring than in the CMHOG 
run. Since the Athena run has a lower resolution at 
i? < 1.1 kpc, on the other hand, it harbors weaker nu- 
clear spirals than in the CMHOG run, as Figure [T] shows. 
The peak density relative to the mean and the pitch angle 
of the spirals a.t R — 0.25 kpc are Spcak/Savg = 3.7 and 
-30° in the Athena run, which can be compared 



with Spcak/S a 



7.9 and ir, 



-34° in the CMHOG 



Figure [2] provides more quantitative comparisons be- 
tween the Athena and CMHOG results. Figure [5^,b 
plot the distributions of surface density and velocities 
along the slit indicated in Figure [TJ The slit starts from 
(x,y) = (0,-1.5 kpc) and ends at (1.9,-2.1) kpc. The 
dust-lane shocks in Athena are slightly weaker and lo- 
cated 0.13 kpc closer to the bar major axis than in 
CMHOG. This is consistent with a smaller nuclear ring 
in the Athena run, since the inner ends of the dust lanes 
are connected to the nuclear ring. The overall veloc- 
ity profiles as well as the amount of the velocity jumps 
from the Athena and CMHOG runs agree well with each 
other. Figure[2j: compares the radial profiles of the mean 
density averaged both azimuthally and temporally over 
t = 300 — 500 Myr. Although the density-weighted ring 
radius is slightly smaller in the Athena run, the difference 
is only ^ 0.05 kpc, again demonstrating good agreement 
between the CMHOG and Athena results. The black ar- 
row marked by Rxi = 4.7 kpc indicates the maximum 
radial extent of the outermost xi -orbit, outside of which 
density is relatively unperturbed. A small spike of the 
mean density around R = Rxi is caused by the trapping 
of inflowing gas from outside. 

Finally, Figure [2]i compares the mass inflow rates for 
bhO and bh7 models as functions of time. Because of 
difference in the way handling fluid variables inside the 



:nctic field configurations (right) of Model bh7MHD01 with A/bh = 
1.5 < a; < —0.5 and y = 2 kpc indicates a slit perpendicular to the 

central hole, the temporal behavior of M at early time is 
quite different between the Athena and CMHOG results, 
although it rapidly converges to a more-or-less constant 
value. In the CMHOG runs, the bhO model has larger M 
at t > 300 Myr than in the bh7 model, while the former 
has slightly smaller M in the Athena runs. Nevertheless, 
M at t > 100 Myr from the Athena and CMHOG runs 
agrees with each other, on average, within a factor of 2, 
implying that the current method of measuring M by 
placing a central hole is reliable. 

To conclude this section, despite the CMHOG and 
Athena codes are quite different in terms of coordinate 
geometry as well as ways of updating fluid variables, the 
results from both codes are in good agreement. This 
confirms not only that the revised force transformation 
in CMHOG used in Paper I is correct, but also that dy- 
namics of a rotating disk can be reliably handled by the 
Cartesian Athena code that we use. 

4. BAR REGIONS 

We now turn to the effects of magnetic fields on bar 
substructures and mass infiow rates. Evolution of the 
gas outside the outermost xi-orbit is quite distinct from 
that inside. In this section we focus on the bar regions 
(i.e., inside the outermost cci-orbit): the changes of mag- 
netic structures in the outer regions will be presented in 
Section [5] 

4.1. Overall Evolution 

Early evolution in the bar regions of our standard 
model 'bh7MHD01 (with /3o = 1 and Mbh = 4x 10^ Mq) 
is not much different from the hydrodynamic counterpart 
in that the bar potential induces crowding of gas orbits 
to create overdense ridges that develop into dust-lane 
shocks. The presence of magnetic fields (1) reduces the 
shock density and (2) provides Maxwell stress for the 
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gas at the shocks, further removing angular momentum 
from it. To illustrate this, Figure|3]plots a surface density 
snapshot together with the magnetic field configuration 
of Model bh7MHD01 at t = 150 Myr. Note the mag- 
netic field lines are approximate proxy to instantaneous 
gas streamlines. Figure |4] plots the variations of gas sur- 
face density, velocities, and magnetic fields along the slit 
at y = 2 kpc shown in Figure [S] Each circle corresponds 
to an individual grid point. At this time, the shocks are 
almost parallel to the y-axis, and are moving toward the 
bar major axis at velocity itgh — 30 km s~^ relative to 
the bar. 

The perpendicular Mach number of the flows along the 
slit relative to the shock is Ai± ~ {ux.i — Ush)/cs — 
13. Here and hereafter, the subscripts 1 and 2 denote 
the preshock and postshock value, respectively. With 
= STTc^Si/i?^ 1 — 8, the usual jump condition for 

one-dimensional stationary MHD shock^ yields the com- 
pression factor r ~ = By^2/By^i = 32.5 and 
Bx,2/Bx^i — 1. Figure m shows that S2/S1 ~ 57 and 
By^2/By^i ~ —40, and i?x,2/Sx,i — 1, indicating that 
the density compression of the dust-lane shocks is more 
than steady one-dimensional shocks allow. More impor- 
tantly, the sign of By is reversed across the shock. This 
is because the shocks are intrinsically two dimensional in 
the sense that the streamlines diverge before the shock 
and bend abruptly to align themselves parallel to the 
dust lanes immediate after the shock. The bending of 
the streamlines at the shock is caused by a strong inflow 
of the gas along the dust lanes from the bar ends. This 
postshock inflow rotates the magnetic fields at the shock, 
making the sign of By reversed. Magnetic tension forces 
from the bent _B-ficlds are highly efficient to remove an- 
gular momentum from the gas moving across shocks, as 
we will show below. 

The gas moving on xi-orbits just inside the outer- 
most xi -orbit collides with that on the 4/1-resonant 
family, subsequently sw it ching to lower a:i-orbits (e.g., 
Contopou los fc Grosb0ll 119801 : lEnglmaier fc Gerhard' 
1997: Paper I). Gas on these lower a:i-orbits is continu- 
ally perturbed by strong total (thermal plus magnetic) 
pressure near the bar ends, and after many revolu- 
tions it eventually transits to a trajectory that crosses 
the dust lane shocks. Figure [5^ plots instantaneous 
streamlines of such gas that passes through Point A at 
(x, y) — (0, 1.5 kpc) before hitting a dust lane in Models 
bh7MHD01, bh7MHD10, and bh7HD at t = 150 Myr. 
The locations of the dust lanes are indicated as thick 
lines. The solid lines in Figure [5b plot the change of the 
rotational velocity along the trajectories shown in 
Figure [5^, while the dotted line draws the rotation curve 
in initial equilibrium. The gas hits the shock at Point 
B and loses a significant amount of angular momentum 
there. Note that in magnetized models experiences 
a dramatic drop at the shock due largely to magnetic 
tension forces from the bent field lines. Coriolis force 
makes the shocked gas rotate faster gradually as it 
moves radially in along the dust lane. In MHD models, 
magnetic stress in the dust lanes continuously removes 
angular momentum from the inflowing gas, causing 

^ For perpendicular isothermal shocks, the density compression 
factor r is a positive root of the quadratic equation r-^ + (1 + /9i )»* — 
PlM?^ = (e.g.. lPriesti[T98^ . 
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Figure 4. Profiles of surface density, velocity, and magnetic fields 
along the slit shown in Figure [S] Note the x- and ^/-components of 
the velocity and magnetic fields are in the direction perpendicular 
and parallel to the shocks, respectively. Each circle corresponds to 
an individual pixel value. 

to increase more slowly and thus making it move further 
in than in the unmagnetized model. 

Figure [5] shows that in models with (3q > 10, the in- 
stantaneous streamlines make an almost closed loop af- 
ter hitting the dust lane at the opposite side, forming 
a nuclear ring in a position where the centrifugal force 
balances the gravity (i.e., where achieves its initial 
value). Due to significant angular momentum loss at 
the shocks, on the other hand, the streamline in Model 
bh7MHD01 with /3o = 1 is closed only after successively 
hitting the dust lane at the opposite side (Point C) and 
then the first one again (Point D), allowing the gas to 
move in very close to the galaxy center. With a large 
supply of gas near the central parts, magnetized models 
would have mass inflow rates greatly enhanced compared 
to unmagnetized models, as will be shown in Section [4.51 

Is the magnetic stress really responsible for further 
infall of the shocked gas toward the galaxy center? 
To check this. Figure [6| plots the radial distribution 
of the angle-averaged, specific angular momentum j = 

(27r)^^ Rv^dcj) relative to the initial value jo as well 
as the the ratio of the angle-averaged Maxwell stress 
Mr,ii, = (Stt^)"^ j^^ BjiBrjydcj) to the mean thermal pres- 
sure P = {2n)-^ J^'' c^^T,d(j) at t = 100 Myr for bh7 
models. Note that j > jo a,t R > 1 kpc for all mod- 
els. This is primarily because gas with high angular 
momentum at large R has moved inward where jo is 
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Figure 5. (a) Instantaneous streamlines of gas that starts 
from Point A {x,y) = (0,1.5 kpc) in Models bh7MHD01 (red), 
bh7MHD10 (blue), and bh7HD (black) at t = 150 Myr. The thick 
lines nearly parallel to the j/-axis represent the dust lane shocks in 
these models. In Model bhTMHDOl, the gas is shocked at Point B, 
subsequently moves in along the dust lane, hits the opposite-side 
dust lane at Point C, and returns to the first one to hit it at Point 
D. (b) Variations of the rotational velocity of the gas along the 
trajectories shown in (a). The initial, equilibrium circular velocity 
is plotted as a dotted line. 

lower If the magnetic stress is important in trans- 
porting angular momentum outward, j would change ac- 
cording to dj/dt ~ —c^Mr^^/P, (e.g.. iBalbus fc HawlevI 
[1998). The peak value of the magnetic stress in Model 
bhTMHDOl is Mr^^/P - 1.7 at i? ~ 1.5 kpc. Assuming 
a constant rate of the angular momentum loss, the net 
change in j expected over an interval of At ~ 100 Myr 
is Aj — 4.4 kpc • km s""'^, which is very close to the dif- 
ference in j between Models bhTMHDOl and bh7HD at 
i? ^ 1.5 kpc. This demonstrates that the magnetic stress 
accounts for further loss of angular momentum and hence 
a smaller ring size in magnetized models. 

The density and magnetic field structures in the bar- 
end regions are quite complicated. Figure [7] shows a few 

^ In our models, a portion of the bulge potential is slowly re- 
placed by the bar potential, with the combined central density 
held fixed. A small mismatch between the shapes of the bar and 
bulge potentials causes the gas in the central regions to temporarily 
expand slightly, leading to j < j'o at i? < 1 kpc. 
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Figure 6. Radial distributions of (a) the azimuthally-averaged 
specific angular momentum j relative to the initial value jo and 
(b) the ratio of the azimuthally-averaged Maxwell stress Mji ^f, to 
the azimuthally-averaged thermal pressure P for all bh7 models at 
t = 100 Myr. 

snapshots of density and magnetic fields in the upper bar- 
end regions with |a:| < 1 kpc and 2 kpc < y < 5 kpc of 
the standard model. The 4/1-spiral shock that forms via 
the collisions of gas on Xi-orbits and on the 4/1-resonant 
family is indicated as the red arrow. Also indicated as 
the white arrow is a dense ridge of gas, called "smudge" , 
created by the con vergence of streamlines at th e backside 
of the dust lanes dPatsis fc Athanassoula|[2000l : Paper I). 
Compressed magnetic fields in the smudge bend abruptly 
at the dust-lane shock, developing a "A" -shaped mag- 
netic wedge that exerts exceedingly large tension force 
{t = 170 Myr). At the same time, the dust lane moves 
toward the bar major axis and becomes shorter in extent. 
A blob of high-density gas at {x,y) ^ (0,3 — 4 kpc) near 
the bar end that was previously supported by the dust- 
lane shock becomes loose at i = 190 Myr, and starts to 
move radially in due to the strong tension force in the 
magnetic wedge. This produces a strong non-steady flow 
that rotates about the galaxy center on xi-orbits just 
outside the dust lanes. This non-steady flow, which is 
stronger in models with smaller /3o, dies out slowly as 
it loses mass near the bar ends when hitting the dust- 
lane or 4/1-spiral shocks, with the lost mass tunneled to 
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Figure 7. Snapshots of logarithm of surface density and the magnetic field configurations near the upper bar-end regions at < 1 kpc 
and 2 kpc < y < 5 kpc of Model bhTMHDOl at (a) t = 170, (b) t = 190, and (c) t = 205 Myr, illustrating the formation of a A-shaped 
magnetic wedge and its detachment from the outermost xi-orbit. In (a), the red and white arrows mark the 4/1-spiral shocks and smudge, 
respectively. 



a nuclear ring along the dust lanes. At about t = 500 
Myr, the overall flow pattern reaches a quasi-steady state 
where fluid quantities do not change much over time. 

4.2. Dust Lanes 

The primary response of gas to the imposed bar poten- 
tial is the formation of dust-lane shocks at the leading 
side of the bar. Of course, the degree of density and in- 
field compression depends on the initial field strength. 
Figure [5] plots the temporal evolution of the peak den- 
sity and magnetic field of the dust lanes measured at 
2/ = 1.5 kpc for bh7 models. Since the dust lanes are 
away from the center, they are unaffected by the BH 
mass. 

The peak density Spoak of the dust lanes starts to rise 
rapidly from t ^ 100 Myr when the bar potential at- 
tains a substantial amplitude. Magnetic fields reduce 
the shock compression, making Speak smaller with de- 



creasing Pq. As the gas moves in along the shocks due to 
angular momentum loss, the bar region inside the out- 
ermost xi-orbit is progressively evacuated, limiting the 
growth of Speak- Models with stronger magnetic fields 
have dust lanes closer to the bar major axis. This makes 
the gas density in models with smaller /3o more centrally 
concentrated, causing Speak to increase for a longer pe- 
riod of time. In the Pq = 1 model, the second increase 
of Speak at t = 200 — 240 Myr is caused mainly by the 
non-steady gas inflows initiated from near the bar ends 
explained above. 

Due to shock compression, dust lanes are also regions 
of strong magnetic fields. Figure [5Jd plots the temporal 
variations of the peak i?-field strength in the dust lanes 
at y ^ 1.5 kpc. In Model bhTMHDOl, Bso^poak reaches 
~ 90/iG before it declines rapidly as the gas and mag- 
netic fields are advected to the central regions. Since 
the dust lanes are extremely narrow with thickness of 
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Table 2 

Properties of Nuclear Rings 



Model Kin (kpc) Ront (kpc) Rnng (kpc) (S)max/So Sring/So (-B3D)max {fJ.G) 

bhTMHDOl 0.03 0.71 0.34 29.5 18.6 14.8 

bh7MHD03 0.24 0.82 0.51 29.4 20.5 20.5 

bhTMHDlO 0.33 0.79 0.56 32.9 21.3 15.1 

bh7HD 057 099 076 34^2 17^5 0.0 

bhOMHDOl 0.08 0.87 0.39 24.7 11.2 26.1 

bh0MHD03 0.25 0.81 0.51 33.8 20.8 25.8 

bhOMHDlO 0.31 0.80 0.56 32.7 21.4 13.4 

bhOHD 0.59 1.07 0.80 32.1 15.2 0.0 



Note. — iJin and -Rout are the inner and outer radii of the ring defined by the positions 
where (E) = (S)niax/5, with (S)max being the maximum density; R^ing is the mass-weighted ring 
radius; Sring is the mean density of the ring; {-B3D)max is the peak strength of magnetic fields 
inside the ring. 
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Figure 8. Temporal evolution of (a) the peak surface density and 
(b) peak field strength in the dust lanes at i; = 1.5 kpc for bh7 
models. Dotted lines in (b) denote the field strength obtained by 
virtual observations with a Gaussian beam width of 0.1 kpc. 

^ 0.02 kpc in our models, the observed field strength of 
dust lanes is likely to depend sensitively on the telescope 
beam size. We illustrate this by plotting as dotted lines 
the field strength after Gaussian smoothing with FWHM 
of 0.1 kpc, which is about ^4 — 7 times weaker than the 
value without smoothing. 

4.3. Nuclear Rings 

Figure ini shows the distributions of surface density and 
magnetic fields in the central 1 kpc regions of all models 
at the end of the runs {t = 800 Myr). Three dotted 
curves in Figure|n^ represent a;2-orbits with x^. = 0.6, 0.8, 
and 1.0 kpc, while those in Figure are for xi-orbits 
with Xc — 0.2, 0.4, 0.6 kpc. It is apparent that except 
for Model bhOMHDOl, the shape of a nuclear ring is well 
described by a;2-orbits. Why is the shape of the ring in 
Model bhOMHDOl at quasi-steady state so different from 
the rest? To address this question, we first present how 



the ring forms and evolves in Model bhTMHDOl, and 
then compare it with the case of Model bhOMHDOl. 

4.3.1. Ring Formation 

Figure [10] plots evolutionary changes in gas morpholo- 
gies and field configurations in the central parts of Mod- 
els bhTMHDOl. As the bar potential grows at early time, 
dust-lane shocks become stronger and tend to move to- 
ward the bar major axis. The amount of angular momen- 
tum loss at the shocks becomes correspondingly larger 
with time. At t = 170 Myr, angular momentum loss at 
the shocks due to magnetic tension is so large that the gas 
moving in along the dust lanes makes a sharp turn near 
the X-axis directly toward the galaxy center. This causes 
the dust lanes to bend into an "L" shape, with the inner 
(outer) ends aligned perpendicular (parallel) to the bar 
major axis. This is unlike in hydrodynamic models where 
angular momentum loss at small R is almost negligible, 
so that most of the gas moving radially in forms a ring at 
R ~ 0.9 kpc. The inner ends of the L-shaped dust lanes, 
which are also shocks, come very close together near the 
center. In Model bh7MHD01, they never make contact 
with each other because of centrifugal barrier provided 
by a central BH. Near the center, gas jumps successively 
from one inner end to the other at the opposite side, and 
thereby keeps losing angular momentum and moving to- 
ward the BH. In a sense, the inner ends of the dust lanes 
act as "a bar within a bar" proposed by iShlosman et al.l 
1 19.891), although the inner bar in our models is a gaseous 
one. 

Both inner and outer dust lanes become stronger with 
time as the bar potential continues to grow until t = 186 
Myr. At this time, angular momentum loss at the outer 
dust lanes is so large that the gas falls in toward the BH 
even before arriving at the x-axis, making the inner ends 
rotate slightly in the clockwise direction (t = 200 Myr). 
The central gas that hits the inner dust lanes is accreted 
to the BH by losing angular momentum further. On the 
other hand, out-lying gas whose orbits do not cross the 
inner dust lanes encircles the center largely by following 
a;2-orbits, gradually transforming to a nuclear ring. As 
the amount of gas accreted to the BH increases, the in- 
ner dust lanes become detached from the ring, forming 
a bar structure within the nuclear ring at i = 300 Myr. 
This inner, perpendicular bar feature persists for ~ 100 
Myr. Since the inner bar is just gaseous structure un- 
supported by any gravitational potential, it decays pro- 
gressively with time, turning into weak trailing spirals 
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Figure 9. Logarithm of the density distribution (color scale) and magnetic field configuration (contours) in the inner 1 kpc regions of all 
models at t = 800 Myr. Dotted curves in (a) draw the a;2-orbits that cut the z-axis at Xc = 0.6, 0.8, 1.0 kpc, while those in (b) are for the 
xi-orbits with Xc = 0.2, 0.4, 0.6 kpc. 

at around t = 400 Myr. These spirals are continuously 
affected by pressure perturbations in the ring, weakening 
slowly with time. 

Similarly to in Model bh7MHD01, the dust lanes in 
Model bhOMHDOl with no BH bend into an "L" shape 
at t = 170 Myr, as shown in Figure [TT] Unlike in the 
former where a central BH provides a strong centrifugal 
force for the inflowing gas, however, the gas in the latter 
is almost unresisted and can flow in directly to the center. 
This allows the inner ends of the dust lanes to merge 
together at the origin {t = 200 Myr). As a large amount 
of gas is accreted to the center during this process, the 
perpendicular parts of the dust lanes start to dissolve 
rapidly with time. 

As in Model bhTMHDGl, the contact points between 



the dust lanes and nuclear ring in Model bhOMHDOl tend 
to move upstream along the outer dust lanes. With the 
inner ends of the dust lanes dissipated away, the gas near 
the central part of Model bhOMHDOl follow highly eccen- 
tric orbits that connect the contact points, resulting in 
an inclined eccentric ring {t — 180 Myr). These eccentric 
orbits are ballistic, implying that gas dynamical effect is 
not significant in the ring. The inclined ring precesses 
slowly in the clockwise direction due to the bar torque 
that tends to align the semimajor axis of the ring paral- 
lel to the bar major axis. The ring in Model bhOMHDOl 
eventually reaches an equilibrium position at t = 480 
Myr where the ring gas follows an a;i-orbits and feels 
zero net torque from the bar. 
Less-strongly magnetized models with /3o > 3 also form 
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Figure 10. Snapshots of logarithm of surface density (left) and 
field eonfigurations (right) in the inner 1 kpc regions of Model 
bh7MHD01 at t = 170, 200, 300, and 400 Myr. 

"L" -shaped dust lanes, with the inner ends connected to 
each other at the center when there is no BH. Because 
angular momentum loss at the dust-lane shocks is not so 
significant as in the /3o = 1 model, however, the contact 
points are always located nearly on the a;-axis. The gas 
orbits connecting the contact points are not so eccentric 
as in the /3o = 1 model, either, with their semimajor axes 
almost parallel to the x-axis. With negligible net bar 
torque, the ring does not precess in these weak B-field 
models. As the ring continuously interacts with the gas 
moving in along the dust lanes, it becomes gradually less 
eccentric and more distributed radially, with its shape 
remaining similar to an X2 -orbit. 

4.3.2. Ring Properties 

Careful inspection of Figure [9] reveals that the ring 
becomes smaller as (3o decreases. The density and mag- 
netic fields in the ring are not uniform, often exhibit- 
ing weak trailing spiral features connected to the lower 
ends of the dust lanes. The pitch angle of these spi- 
rals, notably in Model bhTMHDOl, is less than 10°. To 
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Figure 11. Snapshots of logarithm of surface density (left) and 
field eonfigurations (right) in the inner 1 kpc regions of Model 
bhOMHDOl at t = 170, 200, and 300 Myr. 

quantity the mean properties of the rings. Figure [T^] 
plots the radial profiles of surface density (E), magnetic 
fields (-Bsd), and plasma (/3) averaged both azimuthally 
and temporally over t = 500 — 800 yr for all bh7 mod- 
els. Table [5] gives the inner and outer radii, i?in and 
^out, of the ring, the mass- weighted ring radius i?ring = 

/^°"' (S)i?rfi?/ /^°"* (S)di?, the peak density (S)max, the 

mecLn density 5jj-ing — /^°"'(S)rfi?/(i?out - Rin), and the 
peak magnetic field strength (i33D)max inside the ring in 
each model. Here, Rm and i?out are defined as the radii 
where (S) = (S)inax/5, as in Paper I. The ring becomes 
smaller and more centrally distributed due to larger an- 
gular momentum loss in more strongly magnetized mod- 
els. Small values of i?in in magnetized models imply that 
a large quantity of gas is brought in by magnetic torque, 
readily available for accretion to the central BH. 

In our models, the accretion of magnetic fields to the 
galaxy center is realized by removing magnetic fields in- 
side a central hole together with including the Ohmic 
dissipation at i? < 40 pc. Figure [l2b.c show that the 
strength of B-fields in the ring ((Bsd) ^ 15 — 20^G, or 
equivalently (/3) ~ 1 — 3) in bh7 models does not vary 
sensitively to the initial /3o since the loss of magnetic 
flux via accretion to the galaxy center is larger in mod- 
els with stronger initial fields. While surface density in 
models with ^ 3 is more or less uniformly distributed 
inside i?out, magnetic fields peak at i? ~ i?ring and be- 
come weaker at smaller R. The decline of (Bsd) inward 
of i?iing was caused by the field accretion as well as the 
Ohmic dissipation occurred at t ~ 200 — 350 Myr (see 
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Table 3 

Time-averaged Mass Inflow Rate and Its 
Dispersion 
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Figure 12. Radial distribution of (a) gas surface density, (b) mag- 
netic fleld strength, and (c) plasma /3 averaged both azimuthally 
and temporally over t = 500 — 800 Myr for bh7 models. 

Section nisi. 



4.4. Nuclear Spirals 

Figure IHl shows that ah of the magnetized models do 
not possess well-defined nuclear spirals at the end of the 
runs. Only Model bh7MHD01 has short spiral-like fea- 
tures that are remnants of the dissipating inner ends of 
L-shaped dust lanes. This is in contrast to the hydrody- 
namic results of Paper I that showed that when the gas 
is unmagnetized and cold with c,j = 5 km s~^, bhO mod- 
els have strong leading spirals, while bh7 models allow 
weak (but recognizable) trailing spirals. The absence of 
nuclear spirals in magnetized models is primarily due to 
the fact that rings are smaller in size and more centrally 
concentrated than in hydrodynamic models. The ring 
material is continuously perturbed by the infalling gas 
at the contacted points. The associated thermal pressure 
and magnetic forces inhibit growth and maintenance of 
any coherent structures in the central parts. 

4.5. Mass Inflow Rates 

It has been widely accepted that a galactic bar is an ef- 
ficient means to transport the disk gas all the way to the 
center to fuel a central BH. Paper I showed that this hap- 
pens only when the gas has an effective speed of sound 
Cs > 15 km s^^. In this case, large thermal pressure at 
the contact points spreads out the ring material by per- 
turbing gas orbits away from a;2-orbits, some of which 
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on eccentric orbits flow in directly to the center. On the 
other hand, the gas with Cs < 10 km s"^ forms a narrow 
nuclear ring that prevents further inflows of the gas to 
the center. 

We have seen earlier that magnetic stress in the dust 
lanes removes a significant amount of angular momen- 
tum from the gas in the dust lanes and thus causes it to 
move in close to the galaxy center, potentially increasing 
the mass inflow rates. To quantify the effect of magnetic 
fields on M, Figure [13] plots temporal changes of M for 
all models. The time-averaged value of M and its stan- 
dard deviation AM over t = 200 — 800 Myr are given 
in Table [3l A sudden increase of M at i 180 Myr 
in /3o = 1 models is due to the direct inflows of the gas 
through the inner ends of the L-shaped dust lanes. This 
happens later in models with weaker fields since it takes 
the gas longer to reach the galaxy center due to lower 
angular momentum loss at the shocks. On average, M 
m. (3q — \ models is larger by more than two orders of 
magnitude than in the unmagnetized counterpart. This 
is of course because nuclear rings in the former are more 
centrally concentrated and have a much larger interior 
density than the latter. 

The presence of a central BH in bh7 models circular- 
izes surrounding gas orbits and thus makes the initial 
increase of M less dramatic than in the bhO counterpart 
where infalling gas can plunge directly into the central 
hole. With more-or-less circular orbits, the gas in bh7 
models flows in continuously to the hole, exhibiting sus- 
tained mass inflows. On the other hand, the rapid mass 
inflows in bhO models consume the neighboring gas al- 
most completely, lowering M temporarily before fresh 
gas is supplied from the outer parts. When the central 
regions are filled in, the mass inflows resume. This makes 
M relatively intermittent in bhO models^ 

To quantify the amount of magnetic fields removed (via 
accretion or Ohmic dissipation), we plot in Figure 1141 
temporal changes of magnetic fiux <I>b(<) = /B • rfR 
relative to the initial value for all magnetized models. 
A sudden drop in $b occurs almost at the same epoch 
as the steep rise of M, after which remains almost 

3 Assuming a slow and steady inflow due to magnetic torque, 
I Beck et all ||2005I '| estimated A/ ~ 1 Mq yr~^ , very close to the 
peak value in Model bhOMHDOl, although mass inflows in our 
models occur mostly by gas on non-steady eccentric orbits rather 
than steady near-circular ones. 
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Figure 13. Temporal variations of the mass inflow rates M for (a) models with Mbh = 4 x lO'^ Mq and (b) models with no BH. 
Magnetized models with /3o = 1 have M larger by more than two orders of magnitude than unmagnetized models. In bh7 models, M is 
sustained, while it becomes relatively intermittent in bhO models. 

As explained in Section |3l the outermost xi-orbit di- 
vides a galactic disk essentially into two dynamically dis- 
connected regions. Gas responses inside it are highly dy- 
namic, exhibiting dust-lane shocks and significant angu- 
lar momentum loss leading to the formation of a nuclear 
ring. With no crossing of closed gas orbits, however, the 
evolution of gas outside the outermost xi-orbit does not 
involve shocks and is thus much less dramatic. In unmag- 
netized models, the gas at i? > R^i experiences a weak 
bar torque and moves radially in only slightly, piling up 
at the outermost xi-orbit. 

The presence of magnetic fields makes the gas re- 
sponses in the outer regions completely different from 
that of unmagnetized cases. Figure [T5] shows snapshots 
of surface density (logarithmic colorscale) and magnetic 
fields (contours) in Model bh7MHD01. The boxes ex- 
tend to 10 kpc on either side of the center. The sohd 
oval in each panel draws the outermost xi -orbit, while 
dotted circles in the lower panels with contours indicate 
the CR. At t = 180 Myr, magnetic fields at i? > R^i are 
relatively unperturbed compared to those inside. They 
closely follow the outermost xi -orbit near R = R^i and 
are slightly compressed along the weak trailing spirals 
emerging from the leading side of the bar ends. 

Note that the field lines are distorted near the CR and 
inward of it. This arises by the combined action of the 
bar potential and background shear in such a way that 
small radial gas motions induced by the bar potential cre- 
ate the radial component Bn from the azimuthal com- 
ponent of the initial magnetic fields, which is then 
stretched by background shear to generate new B^j, that 
is going to turn to Bji again by the radial motions. The 
whole process makes a closed cycle, resulting in a secu- 
lar growth of magnetic fields in the outer regions. This 
is an ideal MHD dynamo, as opposed to the mean-field 
dynamo, that does not require the parametric terms in 
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Figure 14. Temporal variations of the magnetic flux for all 
models. The decrease of <l?s relative to the initial values in bh7 
and bhO models is about 17, 8, 1% and 18, 3, 0.3% for /3o = 1, 3, 
10 models, respectively. 

constant. This is because magnetic fields in the central 
0.1 kpc regions remain very weak after the rapid accre- 
tion phase, so that the magnetic flux carried in by the 
gas afterward is not appreciable. The decrease of $b in 
bh7 models amounts to ^ 17, 8, and 1% relative to the 
initial value for models with /3q = 1,3, and 10, respec- 
tively, suggesting that the loss of magnetic flux can be 
significant depending on the field strength. 

5. OUTER REGIONS 
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Figure 15. Snapshots of logarithm of gas surface density (color scale) and magnetic field configurations (contours) in the 10 kpc regions 
of Model bhTMHDOl at t = 180, 300, 450, 600, and 800 Myr. The solid oval in each panel draws the outermost xi-orbit, while the dotted 
circle in the lower panels with contours marks the CR. The rectangular section in (b) is enlarged in Figure [TTl The MHD dynamo occurring 
near the CR produces magnetic arms at t = 450 Myr. Magnetic fields reconnect at the base of the magnetic arms, making density and 
B-fields chaotic in the outer regions at t 600 Myr. 
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Figure 16. Logarithm of (a) gas surface density and (b) magnetic 
field strengtti of Model bh7MHD01 at t = 450 Myr in tlie In 
plane. The azimuthal cut-profiles of E (dotted) and B (solid) at 
-R = 6 kpc are given in (c). 

the induction equation. In our models, the MHD dy- 
namo occurs naturally under a combined action of the 
bar potential and galactic shear. Figure [13] shows that 
the MHD dynamo operates most efRciently near the CR 
regions along the bar minor axis. 

Through magnetic stress, the distorted field lines trans- 
port angular momentum from inside to outside, causing 
gas and newly generated field lines at smaller (larger) R 
to move radially inward (outward). Since the field lines 
do not move across the outermost-xi orbit, those mov- 
ing in are continuously accumulated there and increase 
the field strength at the base of the trailing spiral arms, 
while those moving outward are swept by galaxy rotation 
and added to the outer parts of the arms {t = 300 Myr). 
Although less apparent in Figure [15b. the MHD dynamo 
is also effective near the bar ends just inside the outer- 
most xi-orbit where both Bji and shear are strong: Bfi 
is stretched by the gas flows just outside the dust lanes, 
generating new that moves radially outward and is 
added to the base of the spiral arms. 

Due to the MHD dynamo action occurring near the 
CR, the spiral arms become increasingly more magne- 
tized. The associated strong magnetic pressure expels 
the gas away from the arms, eventually forming "mag- 
netic arms" that are characterized by stronger mag- 
netic fields with lower density than the surrounding re- 
gion at i = 450 Myr. This confirms the results of 
iKulesza-Zvdzik et al.l (|2009l I2O1O0 who first found that 



the bar potential develops magnetic arms in outer re- 
gions. Figure [16] plots gas surface density and mag- 
netic field strength in the In R-cj) plane as well as the 
their cut profiles at i? = 6 kpc for Model bhTMHDOl at 
t = 450 Myr. It is apparent that the magnetic arms at 
R = 4.5 — 7 kpc are approximately logarithmic in shape 
with a pitch angle of ~ 18°. The peak B-field strength 
of the magnetic arms at i? = 6 kpc is Bpoak/^o — 17, oc- 
curring at the location where gas surface density is min- 
imized at Smin/So = 0.1. The gaseous arms produced 
by the expelled gas from the magnetic arms have the 
peak density Spoak/^o — 24 and lead the magnetic arms 
by ~ 15° in the azimuthal angle. The magnetic arms 
keep growing in strength as the MHD dynamo continues 
operating. Since the added fields are preferentially az- 
imuthal, the outer ends of the arms curl back in (t — 600 
Myr) and touch the bar end at the opposite side {t — 800 
Myr). 

Figure [15b shows that surface density and magnetic 
fields in the outer regions of Model bh7MHD01 become 
quite chaotic at the end of the run. This is caused by re- 
connection of magnetic fields occurring at the base of the 
magnetic arms. Figure [TT] illustrates how magnetic fields 
reconnect in the rectangular section shown in Figure [TSb 
from t = 300 to 450 Myr. The shaded (unshaded) regions 
in the upper panels represent domains with positive (neg- 
ative) B^, with the interfaces representing the reversal of 
the field direction. The thick blue line in the upper pan- 
els indicates a part of the outermost xi-orbit. The lower 
panels plot the profiles of B^ along the y = cut. At 
t = 300 Myr, the shaded region at the right side of the 
outermost si-orbit corresponds to a segment of the mag- 
netic arms. The sheared magnetic fields with negative 
Bfj, to the left of the magnetic arms are pushed outward 
as the inner regions expand due to the pressure gradi- 
ent {t = 350 Myr). With the base of the magnetic arms 
acting as a rigid wall, this squeezes the magnetic fields 
of opposite polarity into a very narrow layer, developing 
a thin current sheet at a; ~ 4 kpc {t = 400 Myr). Non- 
zero numerical resistivity allows the fields to reconnect 
in the current sheet through a tearing-mode instability, 
producing numerous magnetic islands distributed along 
the base of the arms {t = 450 Myr). Magnetic islands 
that are also sites of density compression move along the 
magnetic arms in the course of galaxy rotation, and inter- 
act with the surrounding gas to produce chaotic density 
and magnetic structures in the outer regions. 

Figure [18] plots the temporal evolution of the mean 
magnetic energy density = i?^/(87r) relative to the 
initial thermal energy density for bh7 models. The solid 
and dotted lines are for an annulus with i? = 5 — 10 kpc 
and its interior at i? = — 5 kpc, respectively. The 
initial rise of Cm in the interior region at f = 50 — 150 
Myr is due to the formation of dust lanes where i3-fields 
are compressed, while the second rise at t ^ 300 — 350 
is due to the MHD dynamo occurring at the bar ends. 
The MHD dynamo operating near the CR is responsible 
for the growth of Cm in the i? = 5 — 10 kpc annulus for 
t < 500 Myr. The corresponding increase of the magnetic 
energy density is a factor of ^ 2 for Pq = 1 models and ^ 
5 for /3o = 10 models. The ensuing decrease of results 
from the magnetic reconnection that not only reduces 
magnetic energy but also makes the outer disk chaotic. 
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Figure 17. Upper panels: evolution of magnetic field configurations in the rectangular section shown in Figure [Tsl The shaded regions 
represent the domains with positive B^, while those with negative are unshaded. Lower panels: profile of B^ along the y = cut. 
Magnetic fields start to reconnect at t = 400 Myr in a thin current sheet located at i: ~ 4 kpc via tearing-mode instability, and develop 
numerous magnetic islands along the base of magnetic arms at t = 450 Myr. 

1. Comparisons between the results of hydrodynamic 
models using the CMHOG and Athena codes show that 
except for small differences in the ring size and location 
of dust lanes, two results are in good agreement with 
each other. The differences are most likely due to the 
fact that the Athena runs with a uniform Cartesian grid 
have higher spatial resolution at R > 1.1 kpc, so that 
the dust-lane shocks are a bit stronger, resulting in larger 
angular momentum loss than in the CMHOG runs that 
use a non-uniform cylindrical grid. On the other hand, 
the CMHOG runs resolve the central regions much bet- 
ter and thus have stronger nuclear spirals than in the 
Athena runs. The close agreement between the results 
from the difference codes confirms not only that the bar 
forces are correctly implemented in the CMHOG code 
used in Paper I, but also that the Cartesian Athena code 
is reliable in handling dynamics of rotating flows. 

2. Our adopted model for the external gravitational 
potential has an outermost xi -orbit that crosses the x- 
and y-axes at Xc = 3.6 and yc = 4.7 kpc, respectively, 
and has a Jacobi energy Ej = —1.24 x 10^(km s""'^)^, 
relative to which gas responses to the bar potential is 
completely different between inside and outside. Inside 
this orbit (i.e., in the bar regions), there exist families of 
closed xi- and 3:2 -orbit that cross each other, so that gas 
on these orbits collides to produce shocks that eventually 
develop into dust lanes and a nuclear ring. As the gas 
loses angular momentum and moves toward the galaxy 
center, the bar regions becomes progressively emptied, 
decreasing the strength of the dust lanes. Outside this 
orbit, on the other hand, there is no closed orbit and 
the bar forces are weak, resulting in much milder gas 



6. SUMMARY AND DISCUSSION 
6.1. Summary 

We run high-resolution MHD simulations using a mod- 
ified version of the Athena code to study the effects of 
magnetic fields on bar substructures and the mass infiow 
rates in barred galaxies. This work directly extends Pa- 
per I in which we explored the case of unmagnetized disks 
using the CMHOG code. Most previous studies on mag- 
netic fields in barred galaxies employed the mean-field 
dynamo, ignoring the back reaction of magnetic fields 
that evolve passively according to the velocity fields ob- 
tained from hydrodynamic runs (see Section 1 for refer- 
ences). On the other hand, | Kulesza- Zvdzik et al.l (|2009l 
120101 ) and iKulpa-Dybel et a l. (2011) performed three- 
dimensional simulations using ideal MHD models, but 
they focused mainly on the formation of magnetic arms in 
the outer regions and their morphological changes due to 
varying rotation frequency and sound speed. While our 
models are two-dimensional, they have about 10 times 
higher resolution in the in-plane direction than the three- 
dimensional models of IKulesza-Zvdzik et al.l ()2010l ). en- 
abling detailed study of the bar and nuclear regions. 

We consider an infinitesimally-thin, isothermal, ro- 
tating, magnetized gas disk with initially uniform sur- 
face density. The magnetic fields are initially purely 
azimuthal and uniform with strength measured by the 
dimensionless plasma parameter /3q. We fix the sound 
speed to Cs = 5 km s^^ and vary /3o as well as the mass 
Mbh of a BH that controls the rotation curve near the 
galaxy center. The main results of the current paper can 
be summarized as follows. 
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Figure 18. Temporal changes of the ratio of magnetic to thermal 
energy densities in an annulus with i?, = 5 — 10 kpc (solid) and 
in the interior region with R < 5 kpc (dotted) in all magnetized 
models. 

responses than in the bar regions. 

3. Even in the presence of magnetic fields, the bar 
regions produce a pair of dust-lane shocks at the lead- 
ing side of the bar and a nuclear ring, just like in un- 
magnetized models. However, magnetic fields make sev- 
eral quantitative changes in the properties of bar sub- 
structures and mass inflow rates. First, magnetic fields 
compressed in dust lanes tend to reduce the peak den- 
sity of the shocks compared to the unmagnetized coun- 
terpart. Second, the post-shock inflows immediate be- 
hind of the shocks along the dust lanes rotate magnetic 
fields abruptly at the shock fronts. The bent field lines 
exert magnetic tension forces to the gas moving across 
the shocks, removing further angular momentum from 
it. This causes the infalling gas to move closer to the 
galaxy center, forming a nuclear ring that is smaller in 
radius and more centrally distributed than in the hy- 
drodynamic model. Third, small centrally-concentrated 
rings in magnetized models destroy coherent perturba- 
tions in the nuclear regions that would otherwise grow 
into nuclear spirals in magnetized models. Fourth, the 
enhanced density in the nuclear regions increases the 
mass inflow rates greatly. For instance, magnetized mod- 
els with I3q = I have a time-averaged mass inflow rate of 
M ~ (5 — 8) X 10~^ Mq yr~^, which is larger by more 
than two orders of magnitude than that in the unmag- 
netized models. 

4. In models with /3o > 3, a rapid loss of angular mo- 
mentum caused by magnetic stress at dust-lane shocks 
causes them to bend transiently into an "L" shape, with 
the lower ends pointing roughly perpendicular to the bar 
major axis. The lower ends of the dust lanes essentially 
play a role of "bars within bars" since the gas moving 
across them keeps losing angular momentum, moving 
toward the galaxy center. Unsupported by any gravi- 
tational potential, this bar-within-a-bar phase does not 
persist, lasting only for about 100 Myr. As the central 
gas is lost via accretion to a BH, the lower ends of the 



dust lanes dissipate gradually, turning to trailing nuclear 
spirals in models with A^bh = 4 x 10^ M© that decay 
with time. 

5. The shape of a nuclear ring at the end of the run can 
be completely different depending on the parameters. In 
our models, all nuclear rings that form follow an X2-orbit 
closely except for Model bhOMHDOl (with /3o = 1 and 
no BH) in which the ring shape is described rather by an 
zi-orbit elongated along the bar major axis. In the latter 
model, there is no centrifugal barrier associated with a 
central BH, so that the lower ends of the L-shaped dust 
lanes come very close together to merge at the center, 
after which they dissolve rapidly with time. The remain- 
ing gas forms a highly eccentric ring that is inclined with 
respect to the bar minor axis. The inclined ring under 
the bar torque precesses slowly in the clockwise direction 
to align its long axis parallel to the bar major axis. 

6. In hydrodynamic models, the regions outside the 
outermost xi-orbits are almost featureless other than 
possessing a pair of weak trailing spirals that emerge 
from the leading side of the bar ends. In magnetized 
models, however, the outer regions are dynamically ac- 
tive involving an MHD dynamo, magnetic arms, and 
field reconnection. The MHD dynamo occurs near the 
CR and bar-end regions where both the bar forces and 
background shear are strong, so that the bar potential 
induces the radial velocity perturbations, producing Bji 
from B^. Background shear then stretches Bfi to gener- 
ate B^ that subsequently turns into Bn due to the radial 
velocity perturbations, closing the loop of the MHD dy- 
namo cycle. The increase of magnetic energy due to the 
MHD dynamo is about a factor of ~ 2 — 5 with a larger 
value corresponding to stronger initial fields. The ampli- 
fied fields move inward/outward radially, adding to fields 
in the trailing spiral arms as well as their base. As the 
arms become more magnetized, the arm gas is expelled 
by strong magnetic pressure, forming magnetic arms that 
have stronger fields but lower density than the surround- 
ing regions. Magnetic fields with opposite polarity pro- 
duced by the MHD dynamo are compressed into a thin 
layer at the base of the magnetic arms. With non-zero 
numerical resistivity, the fields start to reconnect in the 
layer via conventional tearing-mode instability, produc- 
ing numerous magnetic islands with large density. These 
magnetic islands propagate along the magnetic arms to 
make the outer regions highly chaotic. 

6.2. Discussion 

Our numerical results show that the density and mag- 
netic fields in dust lanes remain strong only for 100 Myr 
around the time when the bar potential achieves the 
full strength. The rather rapid decline of the strength 
of dust lanes is primarily due to the fact that the 
gas located outside the outermost xi-orbit is trapped 
there and unable to move further in. As the gas in- 
side this orbit experiences shocks and falls in toward 
the center, the bar regions become evacuated and the 
density in the dust lanes drops accordingly. Obser- 
vations indicate that most barred galaxies without in- 
ner rings possess a pair of p rominent dust lanes (e.g., 
iKormendv fc Kennicutt)l2 004l th at are nearly stra ight in 
strong bars (|Knapen et al. 2002: [Comeron et al.l [2009) . 
Our results suggest that barred galaxies with strong dust 
lanes should have a stellar bar that is dynamically young. 
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or there should be mechanisms replenishing the gas in the 
bar regions. Candidate mechanisms for the latter include 
a spiral-arm potential outside the CR that perturbs gas 
orbits to send it in across the outermost xi -orbit and in- 
falls of halo gas in the for m of galactic fountains (e.g., 
iFraternali fc Binnev. ,2006, 200l . 

Radio pol arization observations o f barred spiral galaxy 
NGC 1097 (IBeck et al.lll99l 120051) reveal that (1) radio 
ridges roughly coincide with the dust lanes in the bar 
region; (2) the direction of the magnetic fields changes 
strongly in the upstream side, making the polarized emis- 
sion almost absent there; (3) the equipartition strength 
Bp of the regular plus anisotropic random magnetic 
fields is ^ 10 — 20/iG in the nuclear ring; (4) Bp ~ 
7 — 12/xG in the dust lanes; (5) the regular fields pass 
through the nuclear ring, with an inclination angle of 

50°, which appears to deviate from the streamlines of 
the gas flow around the ring; and (6) the regular fields 
in the nuclear region have a spiral shape with large pitch 
angles. Points (l)-(3) are entirely consistent with our 
numerical results. Point (4) appears consistent with our 
results provided that the bar in NGC 1097 is dynami- 
cally young; otherwise density and magnetic fields of the 
dust lanes would not be strong enough to be observed. 
Again, rejuvenation of the dust lanes due to gas inflows 
by spiral arms and/or galactic fountains would be an 
alternative possibility. Points (5) and (6) appear incon- 
sistent with our results that show that magnetic fields 
follow the nuclear rings very closely with a pitch angle 
less than 10° and nearly circular in the central parts (see 
Fig. E]). These discrepancies may owe to the fact that 
our models are limited to a two-dimensional razor-thin 
geometry with purely azimuthal initial fields, and thus 
are unable to capture the potential effects of poloidal 
fields and other dynamical processes that may be im- 
portant in three dimensions. In addition, star formation 
and associated small-scale turbulent dynamo occurring 
in nuclear rings of real galaxies may considerably affect 
the field configurations in the central regions. 

The concept of " bars w ithin bars" was introduced 
by iShlosman efall (|1989( ) to overcome the difficulty, 
caused by the formation of a narrow nuclear ring, of 
a single large-scale bar in bringing the disk gas all 
the way to within ~ 1 pc from the galaxy center to 
feed an AGN. Near- infrared and/or CO observations 
of nearby barred galaxies indeed show that a substan- 
tial fraction of disk g alaxies exhibit double bars (e.g.. 



Shaw et al. 1995; Fri edh et all [ 
1999t iGarcia-Burfllo et all 119981: 

20001 ) ■ with the extent of the secondary (nuclear) bars 
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typically 5 to 7 times shorter than that of the primary 
bar. The origin of secondary bars is yet unclear. They 
may arise from gravitation al instability of a disk contain- 
ing both stars and gas (e.g..lFr edli fc Benzll993HCombesl 



IT99llHeller fc Shlosmanfl994D . or be due to crowding of 
stellar orbits per turbed by gravity of a nuclear ring (e.g., 
iShaw et al]|1993() , or be produ ced by an unstable gas eous 
disk in the central part fe.g.. IShlosman et al1ll989l : see 
also lAthanassoulal [20001 and references therein) . Us- 
ing a dynamically- possible doubly-barred galaxy model, 
IMacieiewski et all (12002,) showed that a hydrodynamic 
response of the gas to a self-consistent, secondary-bar po- 
tential is much weaker than that to the primary bar and 
is thus unable to enhance the mass inflow rates much. 



Our numerical results show that magnetic fields natu- 
rally produce, albeit transiently, a bar-within-bar struc- 
ture that greatly enhances the mass inflow rates even 
without invoking the secondary stellar bar potential. In 
our models, the secondary bar (or inner dust lanes) is 
almost corotating with the outer primary bar and lasts 
only for ^ 100 Myr before turning to weak trailing spi- 
rals. Since we include the effect of stars via a fixed grav- 
itational potential, we are of course unable to consider 
the back reaction of stars to the gaseous gravity. 

Our models show that MHD dynamo action due to 
a bar potential combined with background shear pro- 
duces magnetic arms in the outer regions, consistent 
wit h the results of iKulesza-Zvdzik et al.l (|2009L 12010( 1 
and iKulpa-Dvbel et al.l ( 2011 ) who reported that mag- 
netic arms drift into interarm regions due to a lower 
angular velocity than gaseous spirals. In most exter- 
nal disk galaxies, magnetic field directions based on po- 
larized synchrotron radiation follow optical spiral struc- 
tures fairly well, with stronger total (regular plus turbu- 
lent) fields inside the arms than outside (e.g.. IBeck et al.l 
ll996HFletcher et al.l[20TTl ). However, there are some ex- 
ceptional galaxie s such as IC 342 (e.g.. iKrause . 1993.) and 
NGC 6946 (e.g.. IBeck fc Hoernesill996t ) that are known 
to have strongest magnetic fields at magnetic spiral arms 
that lie in between optical arms. Proposed mecha- 
nisms for magnetic arms include M HD dens ity waves 
in two dimensions (£an fc Lou 199 6: ILou fc Fan 19981) 
and mean-field turbulen t dynamo (|Mos£lll998l : iShukurovl 
[19981 IRohde"etaI1ll999l) . The MHD wave theory invokes 
a systematic phase-shift between fast and slow MHD 
modes res ponsible for optical and magnetic arms, respec- 
tively (L ou fc Faiil Il998). although this requires rigid- 
body rotation over a wide range of radii. In addition, 
MHD waves are subject to a buoyant in stability when th e 
vertical degree of freedom is allowed ([Shukurovl |2005[) . 
On the other hand, the mean-field dynamo theory ar- 
gues that optical arms have stronger turbulent motions 
hence larger turbulent magnetic diffusivity, resulting in 
weaker fields than interarm regions, although it relies on 
uncertain parameters including the dynamo number and 
ignores the dynamical effects of magnetic fields on the 
gas. 

Can the MHD dynamo occurring in the outer regions 
of our numerical models account for observed magnetic 
arms in IC 342 and NGC 6949? We think this is un- 
likely since the connection between bar-induced mag- 
netic arms and those observed is quite uncertain. First 
of all, the disk galaxies with observed w ell-defined mag- 
netic arms are not strongly barred (e.g-. lRegan fc Vogel 
119951 ). In addition, the spiral arms in NGC 6949 have 
weak, m — A modes, whereas the bar potential preferen- 
tially amplifies m = 2 arms. Moreover, it is uncertain if 
the MHD dynamo in three dimensions occurs similarly 
to in two-dimensional in-plane geometry, studied in this 
work. Nevertheless, the MHD dynamo as a mechanism 
for field amplification in barred galaxies is quite attrac- 
tive in that it uses only natural ingredients (bar potential 
and galactic shear) without making any assumption. It 
may have something do with relatively weak magnetic 
arms observed in o uter regions of barred galaxies such 
as NGC 1365 (Be ck et al.ll200l . NGC 1097 and NGC 
1365 (IBeck et al.ll2(105i r although it is questionable if the 
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gaseous arms in these galaxies are driven solely by a bar 
potential. While our models considered only a bar poten- 
tial as a perturbing agent for simplicity, real galaxies also 
have spiral arms in the regions outside a bar. It would 
be interesting to study how the MHD dynamo due to 
a bar potential conspires with a spiral-arm potential to 
generate and shape magnetic fields in outer regions. 
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APPENDIX 

NUMERICAL VISCOSITY AND MAGNETIC DIFFUSIVITY OF THE ATHENA CODE 

In this Appendix we evaluate numerical viscosity and magnetic diffusivity of the Athena code by using the damping 
rates of traveling magnetosonic waves. We begin by writing the MHD equations for isothermal, viscous, and resistive 
plasma 

|j+V-(pv)=0, (Al) 

— +vVv = --^Vp+ (V X B) X B + i^V^v, (A2) 

ot p iirp 

— = V X (v X B) + TyV^B, (A3) 

where v and rj are the coefficient of kinematic viscosity and magnetic diffusivity, respectively. Other symbols have 
their usual meanings. 

For the purpose of measuring the numerical values of v and rj of the Athena code used in this work, we limit ourselves 
to magnetosonic waves propagating along the direction perpendicular to the initial magnetic fields Bq = Boy through 
an initially static, uniform medium with density po. We consider small-amplitude perturbations pi(x, i), vi — vi{x, t)x, 
and Bi — Bi{x,t)y to the density, v eloci ty, and magnetic fields, respectively. Assuming Ipil/po ^ 1, |wi|/cs <JC 1 and 
\Bi\/Bi^ <C 1, we linearize equations (jAip - (jA3|) to obtain a set of the perturbed equations 

dpi dvi 

dvi _ cl dpi Bq dBi ^ ^ d'^vi ^^^^ 
dt pq dx 47r/9o dx dx"^ ' 
dBi dvi d^Bi 

We now seek for the plane- wave solutions pi,vi, Bi cx exp(ifca; — icot) with wavenumber k and frequency w. Plugging 
these into equations (|A4[) - (jA6l) and eliminating pi and Bi in favor of wi, we obtain the dispersion relation 

+ i{v + ri)k^uj'^ - [{cl + v\)k^ + vr]k'^]uj - iriclk^ = 0, (A7) 

where wa = -Bo/(47rpo)^^^ is the Alfven speed. It is trivial to show that when v — rj — Q, equation (|A7|) is reduced to 
the relation 

uj^ ^^l^{cl+vl)k\ (A8) 
for magnetosonic waves. The corresponding (real) eigensolutions are 

sin(fcx — ojoOi (^9) 

where is the amplitude of the perturbations in the absence of viscosity and magnetic diffusivity. 
When u and 77 are non-zero but small (i.e., v,r] <^ cjo/fc^), equation (jA7p yields 

u = uo-i[i^ + vl/{ci + vl)rj]ky2, (AlO) 




implying that viscosity and magnetic diffusivity cause the perturbation amplitude in equation (jA9p to decay exponen- 
tially with time as 



A{t) ~ Ao exp 



(All) 



To evaluate v and 77 separately, we first consider pure sound waves with wavelengt h A = 27r/fc in an unmagnetized 
medium, and initialize in the Athena code the perturbations according to equation (|A9p at t = 0. We fix the initial 
amplitude to Aq = 10"'^, and monitor the temporal decay of the wave amplitude. The resulting decay rates for various 
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Figure 19. Decay rates of acoustic and MHD waves in the Athena test runs as functions oi N = A/Ax, the number of grid points per 
wavelength. 

runs with differing N = X/Ax, the number of grid points per wavelength, are plotted in Figure[Tn]as open diamonds. 
The best fit to the numerical viscosity is found to be 

Note that f„ is proporti onal to Ax^ /X^, c ompletely analogous to the behavior of the numerical conductivity of the 
Athena code reported bv lKim et al.l (|2008[ ). 

Next, we set up magnetosonic waves using equation (|A9|) at t = together with Ao = 10~^ and Cs = va- Figure 
1191 plots as filled circles the numerical damping rates of the magnetosoni c wav es due to both viscosity and magnetic 
diffusivity. After correcting for the viscous contribution using equation (jA12|) . the residual damping rates are fitted 
by the numerical diffusivity as 

- ^ (A) (t^)' (lib)" 

for VA ^ Cs- Note that ry„ and are of the same order. While r/„ is quite small, the related magnetic diffusion 
time Tmag = Ax'^ /r]n can be comparable to the dynamical time (e.g., bar orbital period) if the field strength changes 
substantially across the grid spacing Ax. For instance, Tmag ~ 100 Myr for Ax = 7.3 pc and A ~ 0.5 kpc typical for 
tearing modes of magnetic reconnection occurring in our simulations. 
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